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A modification of the saturation model of deep inelastic scattering at 
small x which includes the Altarelli-Parisi (DGLAP) evolution is presented. 
Significant improvement of the description of the structure function F2 at 
large Q 2 is achieved and a good description of diffractive data is preserved. 

PACS numbers: 13.60.Hb, 12.38.Bx 

1. Introduction 

The saturation model [1] has provided a successful description of HERA 
deep inelastic scattering data. This includes both the F2 data, with the tran- 
sition into the nonperturbative photoproduction region, and the DIS diffrac- 
tive data [2]. Whereas the formulae are particularly appealing through their 
simplicity, they also have an attractive theoretical background: the idea of 
parton saturation. Despite its success, the model suffers from the lack of 
proper scaling violation, i.e. at larger values of Q 2 it does not exactly match 
with QCD evolution described by the DGLAP equations. This becomes 
clearly visible in the region Q 2 > 20 GeV 2 where the model predictions are 
below the data on F2. Therefore, one expects that QCD evolution should 
enhance the cross section in this region. 

We present a modification of the saturation model which preserves its 
success in the low-Q 2 (transition) region, while incorporating QCD evolu- 
tion in the large-Q 2 domain [3] . Since the energy dependence in the large-Q 2 
region is mainly due to the behaviour of the dipole cross section at small 
dipole sizes r, our changes will affect mostly the small-r region. At the 
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same time, we leave the large-r behaviour of the dipole cross section practi- 
cally unchanged which allows to retain a good description of DIS diffractive 
cross section. Recent attempts [4] along the same lines indicate that diffrac- 
tion provides a highly nontrivial restriction on possible modifications of the 
saturation model. 



2. The saturation model and its modification 

We start with a brief review of the saturation model [1]. Within the 
dipole formulation of the j*p scattering, the cross section 

°t,l( x iQ 2 ) = J ^ rdz i>T,L(Q,r,z) a(x,r) tl> T ,L(Q,r,z), (1) 

where *Pt,l are the virtual photon wave functions with the transverse and 
longitudinal polarisation. In the saturation model the following form of the 
dipole cross section a is proposed 

a(x,r) =a {l- exp (-r 2 /AR 2 {x)) } , (2) 

where Ro(x) is the saturation scale which decreases when x — > 0, 

*o(*) = 7^72 (&) X ■ (3) 



GeV^ \x 

In order to be able to study the formal photoproduction limit, the Bjorken 
variable x = xb was modified to be 

/ 4m 2 q \ Q 2 + Am 2 
x = x B 1 + — £\ = — , (4) 



Q 2 J W 2 

where m q is an effective quark mass and W denotes the -y*p center-of-mass 
energy. The parameters of the model, cfq = 23 mb, A = 0.29 and xq = 3-1CT 4 
(for fixed m q = 140 MeV), were found from a fit to small-x data [1]. 

As it is well known [5], in the small-r region the dipole cross section is 
related to the gluon density obeying the DGLAP evolution 

a(x,r) ~ —r 2 a s xg(x,fi 2 ), (5) 

where the scale fi 2 ~ C/r 2 for r — > 0. The equation (5) is valid in the double 
logarithmic approximation in which the constant C is not determined. Ex- 
panding the exponent in eq. (2) for r <g 2Rq(x), we find the gluon density 
in the saturation model 
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For fixed a s this gluon density is clearly scale independent. Thus, we have 
to modify the small-r behaviour of the dipole cross section to include the 
DGLAP evolution and, at the same time, keep the large-r behaviour un- 
changed. This will preserve the idea of saturation, which reflects unitarity, 
and allows a good description of the diffractive cross section. 

Therefore, we propose the following modification of the model (2) 

a(x,r) = a < 1 - exp I — I } , (7) 

where the scale fi 2 is assumed to have the form 

» 2 = C/r 2 + »l (8) 

The parameters C and /x 2 , wm be determined from a fit to DIS data. The 
scale hq S> A allows to freeze the value of the gluon distribution for large r 
at a perturbative scale which leads to &(x, r) ~ do, as in the original model. 
The transition between small and large r depends on x but in detail it might 
be different from the original formulation. Thus, the modified form mimics 
in a consistent way the saturated behaviour of the dipole cross section. 

In a first approximation, g(x,/j, 2 ) is evolved with the leading order 
DGLAP equation in which quarks are neglected in the spirit of the small- x 
limit. The starting gluon distribution at the initial scale Qq = 1 GeV 2 

xg{x,Q 2 ) = A g x-^(l-xf-\ (9) 

where A g and X g are another fit parameters. The exponent determining the 
large x behaviour is motivated by the recent MRST parameterisation [6]. 

For small r, the exponential in (7) can be expanded in powers of its 
argument, and relation (5) with the running a s = a s (/j, 2 ) is found. In con- 
trast to the model (2), the rise in 1/x now has become r-dependent. When 
inserting a into (1) and convoluting with the photon wave function, the 
integrand peaks near r ~ 2/Q for large Q 2 , and the argument of the gluon 
density turns into n 2 ~ Q 2 . Consequently, with increasing Q 2 , DGLAP 
evolution will strengthen the rise in 1/x, whereas in the original saturation 
model the power of 1/x had been constant 

3. Fit results and comparison with data 

We performed a global fit to the DIS data with x < 0.01 in the range 
0.1 < Q 2 < 500 GeV 2 . For the HERA experiments the new 1996-97 data 
from HI and ZEUS were used [7] together with the E665 experiment data 
[8] (in total 330 points) The statistical and systematic errors were added 
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Fig. 1. The dipole cross section as a function of the dipole size r for different values 
of log 10 x = —2,— 3, ...,—7 (from the left to the right). The original saturation 
model: dashed lines and the improved model: solid lines. 

in quadrature. For a full discussion of fit details see [3]. With the fixed 
quark mass m q = 140 MeV, the value of X 2 l^df = 1-18 was found (for 
the original model refitted to the new data X 2 l^df — 3) with the following 
fit parameters: C = 0.26, //§ = 0.52 GeV 2 , A g = L2 and X g = 0.28. In 
addition, for the agreement with the diffractive data we fix the normalization 
of the dipole cross section to the original saturation model value do = 23 mb. 

The form of the new dipole cross section (7) is shown in Figure 1 (solid 
lines) for different values of x. As expected the main modification in com- 
parison to the model (2) (dashed lines) lies in the small-r region. In Figure 2 
the global characteristic of the data description is shown. Namely, we plot 
the effective slope A(Q 2 ), obtained from the parameterisation of F2 at small 
x, F2 ~ x~ A (<2 ) ) for fixed Q 2 . For the shown curves, A was computed using 
the relation F2 = Q 2 / {A-K 2 a em ) o^-fL with the two discussed forms of dipole 
cross sections. As expected, the inclusion of the DGLAP evolution for small 
r significantly improves agreement with the data at large Q 2 while at small 
Q 2 the results are practically the same. 

An important aspect of the dipole models is their straightforward de- 
scription of diffractive processes. In particular, the constant ratio of the 
inclusive over diffractive cross sections as a function of x finds a natural ex- 
planation in the saturation model [1, 2]. In DIS diffraction the cross section 
is dominated by the contribution from large dipole sizes r. Since the large-r 
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part of the dipole cross section is practically unchanged in our modification, 
the description of diffractive data is as good as in the original saturation 
model. The only change introduced is a different treatment of the colour 
factors for the qqg component, see [3] for more details. 
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